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In order to evade the proton decay which appears when the Zee model is embedded into a SUSY GUT scenario 
with _R-parity violation, a new idea based on a discrete Z2 symmetry is proposed. Under the symmetry Z2, the 
quark and lepton mass matrices are tightly constrained. The admissible form of the radiatively-induced neutrino 
mass matrix is investigated. 



1. Introduction 

The Zee model Q] is one of promising models 
of neutrino mass generation mechanism, because 
the model has only 3 free parameters and it can 
naturally lead to a large neutrino mixing g], es- 
pecially, to a bimaximal mixing However, the 
original Zee model is not on a framework of a 
grand unification theory (GUT), and moreover, 
it is recently pointed out Q that the predicted 
value of sin 2 28 so i ar must be satisfied the relation 
sin 2 29 so i ar > 0.99 for Am 2 solar / Am 2 atm ~ 1CT 2 , 
i.e., 

The conclusion cannot be loosened even if we take 
the renormalization group equation (RGE) effects 
into consideration. 

The simple ways to evade the constraint (1.1) 
may be as follows: One is to consider || that 
the Yukawa vertices of the charged leptons can 
couple to both scalars <pi and <fi 2 - Another one 
Js) is to introduce a single right-handed neutrino 
vr and a second singlet Zee scalar S + . Also, a 
model with a new doubly charged scalar k ++ is 
interesting because the two loop effects in such 
a model can give non-negligible contributions to 

*E-mail address: koide@u-shizuoka-ken.ac.jp 

^E-mail address : ambar@anp.saha.ernet.in 

* Present address: Saha Institute of Nuclear Physics, 1/AF 

Bidhannagar, Kolkata 700 064, India. 



the neutrino masses M. As another attractive 
model, there is an idealq] that in an i?-parity vio- 
lating supersymmetric (SUSY) model we identify 
the Zee scalar h + as the slepton e#. Then, we can 
obtain additional contributions from the down- 
quark loop diagrams to the neutrino masses, so 
that such a model can be free from the constraint 
(1.1). 

However, these models have not been embed- 
ded into a GUT scenario. As an extended Zee 
model based on a GUT scenario, there is, for 
example, the Haba-Matsuda-Tanimoto model ||. 
They have regarded the Zee scalar h + as a mem- 
ber of the messenger field M±q + M\q of SUSY- 
breaking on the basis of an SU(5) SUSY GUT. 
However, their model cannot escape from the con- 
straint (1.1) because the radiative masses are only 
induced by the charged lepton loop diagrams. 

In the present paper, we will investigate an 
extended Zee model which is based on a frame- 
work of a SUSY GUT with i?-parity violation, 
and which is free from the severe constraint (1.1). 
Usually, it is accepted that SUSY models with Im- 
parity violation are incompatible with a GUT sce- 
nario, because the i?-parity violating interactions 
induce the proton decay. In order to suppress the 
proton decay due to the i?-parity violating terms, 
we will introduce a discrete Z2 symmetry. 

2. How to evade the proton decay 

We identify the Zee scalar h + as the slepton 
which is a member of SU(5) 10-plet sfermions 
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tpio . Then, the Zee interactions correspond to the 
following i?-parity violating interactions 

A*#f)f(%)f(<M*A S 

-[(ei)i(^)j - ("IW^jKe^k 
+ W L Ud R )^ ~ {d R )^{v L ) ]{d L ) ka ) , (2.1) 

T 

where ip c = Cip and the indices ■■■), 
(A, B, ■ ■ •) and (a, (3, ■ ■ •) are family-, SU(5)gc/t- 
and SU(3) co / or -indices, respectively. However, in 
GUT models, if the interactions (2.1) exist, the 
following i?-parity violating interactions will also 
exist ' 

A*#f)f(V>ioW^)f 

-MWL)j-^ L W L )j](e%) k 
-[(elU^)^ (d R )t(e L Uu L ) ka 

+ WlUA)1 - (dR)TPL)j]{dL)ka} , (2.2) 

which contribute tojhe proton decay through the 
intermediate state d R . 

In order to forbid the contribution of the in- 
teractions (2.2) to the proton decay, for example, 
we can assume that the i?-parity violating inter- 
actions occur only when the field ipi of the third 
family is related, i.e., we assume the interactions 

>?iM)t WioWW , (2.3) 

instead of the interaction (2.2). Then, the terms 
^ij(d R )i(d\)j {u R )3 cannot contribute to the pro- 
ton decay. In order to realize the constraints 

Ar2 = 4s = X k 31 = for k = 1,2 , (2.4) 

we introduce a discrete symmetry Z 2 , which ex- 
actly holds at every energy scale, as follows: 

(V>g)i -» , (^)i -> Vi(^k)i , 

(V'io)i -» &(Vuo)i , (V'io)i -» ^(V'io)i , 

(2.5) 



where r\i and £j take 

r?=(+l,+l,+l) , C = (-1,-1,+1) , (2.6) 

under the Z2 symmetry. Then, the Z2 invariance 
leads to the constraints (2.4). 

However, if the RGE effects cause a mixing be- 
tween the first and third families, the interactions 
(2.3) can again contribute to the proton decay. If 
we assume that 5 and 5 Higgs fields H u and Hd 
transform as 

H u — > +H U , Hd — > +Hd , (2.7) 

under the Z 2 symmetry, the up-quark mass ma- 
trix M u is given by the form 

(c u d u \ 
d u b u . (2.8) 
a u J 

This guarantees that the top quark u 3 in the R- 
parity violating terms (2.3) docs not mix with 
the other components {u\ and 112) even if we take 
the RGE effects into consideration, so that the 
interactions (2.3) cannot contribute to the proton 
decay at any energy scales. 

On the other hand, the down-quark mass ma- 
trix Md and the charged lepton mass matrix M e , 
which are generated by the Higgs scalar Hd, have 
the form 

/ \ 
M d = Mj = 1. (2.9) 
V a i a 2 a 3 J 

The mass matrix form (2.9) cannot explain the 
observed masses and mixings. In order to give 
reasonable masses and mixings of the quarks 
and charged leptons, we must consider additional 
SU(5) 45-plet Higgs scalars, which do not con- 
tribute to the up-quark mass matrix because 
tp 10 M u tpf belongs to (10 x 10) sym metric- Then, 
we obtain the down-fermion mass matrices 

/ c[ c> 2 c' 3 \ 

M d = b\ b' 2 b' 3 , (2.10) 

\ ai + a[ a 2 + a' 2 a 3 + a 3 J 

( -3c'! -3c 2 -3c' 3 \ 
M T e = -3b[ -3b' 2 ~3b' 3 

\ ai — 3a[ a 2 — 3a 2 a 3 — 3a 3 J 

(2.11) 
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where a! i and cj) denote contributions from 
the 45-plet Higgs scalars and H^/ which 

transform H { ^ ] -► +H { + ] and -> -Hi~ ) 

under the symmetry Z2, respectively. Note that 
if we consider either or we cannot still 

give realistic down-fermion masses. For example, 
if we have H$ (H u and H d ) and H^ 5 \ we obtain 
an unwelcome relation 

2 1 2 1 2 / 2 1 2 1 2^ 
m r + + m e — (m b + m s + m d ) 

-8^(|6:| 2 + |c:| 2 )>0, (2.12) 

from the trace of M\M e - M d M\. Therefore, in 
order to obtain realistic down-fermion mass ma- 
trices, we need, at least, the three types of the 
Higgs scalars H 5 , H^' and . 

However, such additional Higgs scalars 
cause another problems. One is a problem of the 
flavor changing neutral currents (FCNC). This 
problem is a common subject to overcome not 
only in the present model but also in most GUT 
models. The conventional mass matrix models 
based on GUT scenario cannot give realistic mass 
matrices without assuming Higgs scalars more 
than two. For this problem, we optimistically 
consider that only one component of the linear 
combinations among those Higgs scalars survives 
at the low energy scale \i — Al (Al is the elec- 
troweak energy scale), while other components 
are decoupled at /i < (Ax is a unification 
scale). 

Another problem is that the 45 Higgs scalars 
can have vacuums expectation values (VEV) at 
the electroweak energy scale A^, so that the Z2 
symmetry is broken at /j, = A^. Therefore, the 
proton decay may occur through higher order 
Fcynman diagrams. In the conventional GUT 
models, it is still a current topic whether the col- 
ored components of the SU(5) 5-plet Higgs scalar 
can become sufficiently heavy or not to suppress 
the proton decay. We again optimistically assume 
that the colored components of the 45-plet Higgs 
scalars are sufficiently heavy to suppress the pro- 
ton decay, i.e., that such effects will be suppressed 
by a factor (A L /A X ) 2 - 
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Figure 1. Radiatively induced neutrino mass 
through the down-quark loop. The vertices A, B, 
C and D are given by (E/fA£/f)y (Ufisn, 
(U^M d U d R ) ku (U$\Ut) ml (Ut) 3k , and 
(UjmlUfyki, respectively. 



3. Radiatively induced neutrino masses 

We define fields Ui, di and as those corre- 
sponding to mass eigenstates, i.e., 

H mass = u L U^M u U%u R + d L U c ?M d U d R d R 

+e L U^M e U e R e R + h.c. , (3.1) 

and fields vli as partners of the mass eigenstates 
e^i, i.e., l^i = (vu^Li)- We define the neutrino 
mass matrix M v as 

H v mass — V^MyVL ■ (3.2) 

Therefore, a unitary matrix U£ which is defined 
by 

Ul T M v Ul = D V = diag(r<, ml) , (3.3) 

is identified as the Maki-Nakagawa-Sakata- 
Pontecorvo |lC| neutrino mixing matrix 
Umnsp = U£. 

In addition to the i?-parity violating terms 
(2.1) and (2.2) [(2.3)], we assume SUSY break- 
ing terms ip^ip w H^ (and ip^wH-^). For the mo- 
ment, we do not consider e c R -H d mixing as in 
the original Zee model. (The H d contribution to 
the neutrino mass matrix will be discussed in the 
end of this section.) Then, the neutrino masses 
are radiatively generated. In Fig. [y, we illustrate 
the Feynman diagram for the case with the down- 
quark loop. The amplitude is proportional to the 
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coefficient 

{U d ^\Ul) l3 {U d L ) Zn -{U^M d U d R ) kl 

\U%\ T Ui) mi {Uih k -{U?mlU d R ) nm 

= (m 2 d XUi) 3l (M d \U e L ) 3j , (3.4) 

where (rh d )ij are coefficients of (d\)i(dft)j , and 
(X)ij — Xfj. Similarly, we obtain the contribu- 
tions from the charged lepton loops. Therefore, 
the radiatively induced neutrino mass matrix M v 
is given by the following form 

(M v )ij = {ftg e 3 + f]gt)K e + (ffgf + ffgf)K d , 

(3.5) 

where Kf (/ = e, d) are common factors indepen- 
dently of the families, and 

ft = (M e T AC/£) 3l , gf = (ml T XUl) 3i , 
ff = (M d XU e L ) 3l , gf = (m 2 XU L ) 3l . 

(3.6) 

On the other hand, the contributions due to 
the Hf - e + mixing are as follows. There are no 
contributions from bilinear terms H u (5)[H d (5) + 
SiV'(5)], because we can always eliminate the 
contributions from H u (5)ip(5) by re-definition of 
the scalar H d (5). Also, there are no contribu- 
tions from H d (5)H d (5)il}(10), because is 
anti-symmetric in SU(5) indices. However, the 
terms H d (h)H d {A5)il) {10) can cause the - e + 
mixing. These contributions affect only to the 
charged lepton loop diagrams. 

The final result which includes the contri- 
butions is given by 

Mij = {ftg'l+fi9l)K tt +<jtg'}+ffgl)K d +FZ j K' e , 

(3.7) 

where 

1% = (UfXMeM'JuDv + (<<-> j), (3.8) 

and M e denotes the contributions due to 
E^e L H+[H+ denotes (H d (E))+, (^ +) (45))+ 

and (Hj~^ (45)) + ]. We have already assumed that 
only one component of the linear combinations 
among those Higgs scalars survives at the low en- 
ergy scale. Then, we can regard M e as M e = M e , 
so that we obtain Fh as follows: 

F% - (X'D 2 e ) ij+ (i ~ j) = X'^mf-mf), (3.9) 



where A = J7| T A?7£ is antisymmetric tensor as 
well as A , and ml = (m e , m M , m T ). Note that, in 
this case, the third term in (3.7) has the same 
mass matrix form as that in the original Zee 
model. 

4. Phenomenology 

In the SUSY GUT scenario, there are many 
origins of the neutrino mass generations. For ex- 
ample, the sneutrinos vn, can have the VEV, and 
thereby, the neutrinos vn acquire their masses 
(for example, see Ref.fLl|). Although we cannot 
rule out a possibility that the observed neutrino 
masses can be understood from such compound 
origins, in the present paper, we do not take such 
the point of view, because the observed neutrino 
masses and mixings appear to be rather simple 
and characteristic. We simply assume that the 
radiative masses are only dominated even if there 
are other origins of the neutrino mass generations. 

However, even under such the assumption, we 
still have many parameters as shown in the ex- 
pression (3.7). For simplicity, we assume that the 
contributions from the charged lepton loop are 
dominated compared with those from the down- 
quark loop, i.e., K e ^> K d , which corresponds to 
the case m 2 (c?L) — m 2 {d R ) 3> m 2 (ei) — m 2 (e^). 
(We can give a similar discussion for the case 
K d ^> K e . In the present paper, we do not discuss 
which case is reasonable.) Furthermore, we ne- 
glect the third term (Ht contributions) in (3.7). 
Then, the neutrino mass matrix (3.5) becomes a 
simple form 

(M„)ij — m (figj + f 3 gi). (4.1) 

Hereafter, for convenience, we will normalize fi 
and gi as 

l/i| 2 + l/ 2 | 2 + |/ 3 | 2 = i, M 2 + bl 2 + bl 2 = i. 

(4.2) 

In the most SUSY models, it is taken that the 
form of rh 2 (/ = e, d) is proportional to the 
fermion mass matrix Mr. Then, the coefficients 
gi are proportional to /,-, so that the mass matrix 
(4.1) becomes {M v )ij — 2mofifj, which is a rank 
one matrix. Therefore, we rule out the case with 
nii oc Mf. 
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For convenience, hereafter, we assume that /j 
and gi (i = 1,2,3) are real. The mass eigenvalues 
and mixing matrix elements for the neutrino mass 
matrix (4.1) are given as follows : 



in 



= (1 + e)77i , 



ml = , 



(4.3) 



U a = 

U a = 

u i3 = 



1 fi 



V2vm ' 
1 fi - gi 

fj9k 



(4.4) 



Sijk- 



where 

£ = fm + $292 + f393 • (4.5) 

As seen in (4.3), the mass level pattern of the 
present model shows the inverse hierarchy as well 
as that of the Zee model. From (4.3), we obtain 



Ara 2 



(m£) 2 - (m^) 2 = -Aeml , 
(m» 3 ) 2 -(m» 2 ) 2 = -(l-e) 2 m 2 , 



R 



4c 



(4.6) 
(4.7) 



^m 2 2 (1 - e) 2 

For a small R, the mixing parameters 
sin 2 28 so i ar , sin 2 29 atm and C/ 2 3 are given by 

sin 2 26 solar = 4U 2 U U 2 2 = j^Ul - Si? , 

(4.8) 

sin 2 29 atm = 4Z7 2 2 3 f/ 3 2 3 



(l- £ 2 ) 



2 , 2 [/2/ 3 + <72<73-e(/352 + / 2 .93)] 2 , (4.9) 



^13 = 1 - 



fi + 9\ - I? f i9i 



(4.10) 



1 -e 2 

Let us demonstrate that the mass matrix form 
(4.1) has reasonable parameters for the observed 
neutrino data. The atmospheric neutrino data 
|l2[ suggests a large v^-Vr mixing, i.e., sin 2 2# 23 — 
1. It is known that the nearly bimaximal mixing 
is derived from the neutrino mass matrix M v with 
2^3 symmetry [O. Therefore, we take 



/l - s a, fi — fz — ~^ c a 7 
91 = C/3, 92 = 93 = -^ S P 



(4.11) 



where c a = cos a, s a — sin a and so on. Then, 
the parameterization (4.11) gives 



e = sin(a — 0) , 

sin 2 29 so i ar = cos 2 (a + (3) , 
sin 2 29 atm = 1 , 



U 







(4.12) 

(4.13) 
(4.14) 
(4.15) 



We assume that the values of a and (3 are highly 
close each other, i.e., sin(a — (3) ~ 10~ 2 . The 
result (4.13) with a ~ j3 is free from the con- 
straint (1.1) in the original Zee model, so that we 
can fit the value of sin 2 29 so iar with the observed 
value HH sin 2 29 so iar ~ 0.8 from the solar neu- 
trino data by adjusting the parameter a (~ /?). 
Of course, the parameterization (4.11) is taken 
only from the phenomenological point of view, so 
that the results (4. 13)- (4. 15) are not theoretical 
consequences in the present model. 

From the recent atmospheric and solar neu- 
trino data 00 R ~ (4.5 x 1CT 5 eV 2 )/(2.5 x 
1CT 3 ) eV 2 ) = 1.8 x 1CT 2 , we estimate s = 4.5 x 
lO" 3 , and m =i m v x ~ |m£| ~ t±m\ tm = 
0.050 eV. The effective neutrino mass (m„) 
from the neutrinoless double beta decay exper- 
iment is given by (m v ) = (M„)n — 2moc a S0 ~ 
movT — sin 2 28 so iar — 2.2 x 10~ 3 eV, where we 
have used the observed value Q| sin 2 29 so iar — 
0.8. 

5. Conclusion 

In conclusion, we have proposed a neutrino 
mass matrix model based on a SUSY GUT model 
where only top quark takes i?-parity violating in- 
teractions and the Z 2 symmetry plays an essential 
role, so that the proton decay due to the i?-parity 
interactions can be evaded safely. 

The general form of the radiatively induced 
neutrino mass matrix in the present model is 
given by the expression (3.7). Since the form 
(3.7) has too many parameters, we have investi- 
gated phenomenology for a simplified case (4.1). 
Although we have investigated the case e ~ [e 
is defined by (4.5)], the case 1 — e ~ is also pos- 
sible. However, the purpose of the present paper 
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is to give the general form of the radiatively in- 
duced neutrino mass matrix (3.7) under the dis- 
crete symmetry Z2. A systematic study of the 
mass matrix (3.7) will be given elsewhere. 
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